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Abstract 

This is the first in a series of two papers in which we derive a T-expansion for 
a two-dimensional non-local Ginzburg-Landau energy with Coulomb repulsion, also 
known as the Ohta-Kawasaki model in connection with diblock copolymer systems. In 
that model, two phases appear, which interact via a nonlocal Coulomb type energy. 
We focus on the regime where one of the phases has very small volume fraction, thus 
creating small "droplets" of the minority phase in a "sea" of the majority phase. In 
this paper we show that an appropriate setting for T-convergence in the considered 
parameter regime is via weak convergence of the suitably normalized charge density in 
the sense of measures. We prove that, after a suitable rescaling, the Ohta-Kawasaki 
energy functional T-converges to a quadratic energy functional of the limit charge 
density generated by the screened Coulomb kernel. A consequence of our results is 
that minimizers (or almost minimizers) of the energy have droplets which are almost 
all asymptotically round, have the same radius and are uniformly distributed in the 
domain. The proof relies mainly on the analysis of the sharp interface version of 
the energy, with the connection to the original diffuse interface model obtained via 
matching upper and lower bounds for the energy. We thus also obtain a characterization 
of the limit charge density for the energy minimizers in the diffuse interface model. 



1 Introduction 

In the studies of energy-driven pattern formation, one often encounters variational problems 



with competing terms operating on different spatial scales [26[|27[|32[|37 ,47 50 52 . Despite 
the fundamental importance of these problems to a multitude of physical systems, their 
detailed mathematical studies are fairly recent (see e.g. [8}{l2}|20|[29}|45j). To a great extent 
this fact is related to the emerging multiscale structure of the energy minimizing patterns 
and the associated difficulty of their description [^ [lT|[T7][30|[33] . In particular, the popular 
approach of T-convergence [5] is rendered difficult due to the emergence of more than two 
well-separated spatial scales in suitable asymptotic limits (see e.g. (9|- |T2^[T7|[30|[33|[46] ) . 

These issues can be readily seen in the case of the Ohta-Kawasaki model, a canonical 
mathematical model in the studies of energy-driven pattern forming systems. This model, 
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originally proposed in 40 to describe different morphologies observed in diblock copolymer 



melts (see e.g. |4j) is defined (up to a choice of scales) by the energy functional 

£ M = I ( £ ^\Vu\ 2 + W(u))dx+l [ [ (u(x)-u)G (x,y)(u(y)-u)dxdy, (1.1) 



where f2 is the domain occupied by the material, u : £1 — > M is the scalar order parameter, 
W(u) is a symmetric double-well potential with minima at u = ±1, such as the usual 
Ginzburg-Landau potential W(u) = j(l — u 2 ) 2 , e > is a parameter characterizing in- 
terfacial thickness, u G (—1,1) is the background charge density, and Go is the Neumann 
Green's function of the Laplacian, i.e., Go solves 

-AG (x,y) = 5(x-y)-^- I G (x,y)dx = 0, (1.2) 

II J r2 

where A is the Laplacian in x and 8(x) is the Dirac delta- function, with Neumann boundary 
conditions. Note that u is also assumed to satisfy the "charge neutrality" condition 
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udx = u. (1-3) 



n 



Let us point out that in addition to a number of polymer systems [16 , 39 49] , this model 
is also applicable to many other physical systems due to the fundamental nature of the 
Coulombic non- local term in ([Tl]) [7][l8}[23](32]j35][38] . Because of this Coulomb interaction, 
we also like to think of u as a density of "charge" . 

The Ohta-Kawasaki functional admits the following "sharp-interface" version: 

E[u] = £ - |Vw| dx + \ j (u(x) - u)G(x, y){u(y) - u) dx dy, (1.4) 
z Jn z Jn Jn 

where G(x, y) is the screened Green's function of the Laplacian, i.e., it solves the Neumann 
problem for the equation 

-AG + K 2 G = 5(x-y), (1.5) 
where k := l/yW"(l) > 0. Note also that in contrast to the diffuse interface energy in 



(1.1), for the sharp interface energy in (1.4) the charge neutrality constraint in (1.3) is no 
longer imposed. This is due to the fact that in a minimizer of the diffuse interface energy, 
the charge of the minority phase is expected to partially redistribute into the majority 
phase to ensure screening of the induced non-local field (see a more detailed discussion in 
the following section). 



The two terms in the energy (1.4) are competing: the second term favors u to be 
constant equal to its average u, but since u is valued in {+1, —1} this means in effect that 
it is advantageous for u to oscillate rapidly between the two phases u = +1 and u = — 1; 
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Figure 1: Two-dimensional multi-droplet patterns in systems with Coulombic repulsion: 
a local minimizer of the Ohta-Kawasaki energy on a rectangle with periodic boundary 
conditions; a local minimizer of the sum of two-point Coulombic potentials on a disk with 



Neumann boundary conditions. Taken from 35 , 42 



the first term penalizes the perimeter of the interface between the two phases, and thus 
opposes too much spreading and oscillation. The competition between these two selects 
a length scale, which is a function of e. In the diffuse interface version (1.1), the sharp 
transitions between {u = +1} and {u = —1} are replaced by smooth transitions at the 
scale e > as soon as £ < 1. 

In one space dimension and in the particular case u = (symmetric phases) the behavior 



of the energy can be understood from the work of Miiller 1 33 1 : the minimizer u is periodic 
and alternates between u = +1 and u = —1 at scale e 1 / 3 (for other one-dimensional results, 
see also |43||55| ). In higher dimensions the patterns of minimizers are much more complex 
and are not well understood. The behavior depends on the volume fraction between the 
phases, i.e. on the constant u chosen, and also on the dimension. When u < 0, we call 
it = —1 the majority phase and u = +1 the minority phase, and conversely when u > 0. In 
two dimensions, numerical simulations lead to expecting round "droplets" of the minority 
phase surrounded by a "sea" of the majority phase (see Fig. [TJ for sufficient asymmetries 
between the majority and the minority phases (i.e., for u sufficiently far away from zero) 
34 35 40 42 . The situation is less clear for u close to zero, although it is commonly 



believed that in this case the minimizers are one-dimensional stripe patterns 13,34,35,40 



In all cases, minimizers are intuitively expected to be periodic. However, at the moment 
this seems to be very difficult to prove. The only general result in that direction to date 
is that of Alberti, Choksi and Otto [I], which proves that the energy of minimizers of the 
sharp interface energy from ( 1.4 ) with no screening (with k = and the neutrality condition 
from (1.3)) is uniformly distributed in the limit where the size of the domain O goes to 



3 



infinity (see also [8j|48]). Their results, however, do not provide any further information 
about the structure of the energy-minimizing patterns. Note in passing that the question 
of proving any periodicity of minimizers for multi-dimensional energies is unsolved even 
for systems of point particles forming simple crystals (see e.g. |30[|46| ), with a notable 
exception of certain two-dimensional particle systems with short-range interactions which 
somehow reduce to packing problems (41 51 53 . Naturally, the situation can be expected 



to be more complicated for pattern forming systems in which the constitutive elements are 
"soft" objects, such as, e.g., droplets of the minority phase in the matrix of the majority 
phase in the Ohta-Kawasaki model. 

Here we are going to focus on the two-dimensional case and the situation where one 
phase is in strong majority with respect to the other, which is imposed by taking u very close 
to —1 as e — > 0. Thus we can expect a distribution of small droplets of u = +1 surrounded 
by a sea of u = — 1. In this regime, Choksi and Peletier analyzed the asymptotic properties 



of a suitably rescaled version of the sharp interface energy (1.4) with no screening in |14 



as well as ( 1.1 ) in 15 1 . They work in the setting of a fixed domain Q, and in a regime where 



the number of droplets remains finite as e — > 0. They showed that the energy minimizing 
patterns concentrate to a finite number of point masses, whose magnitudes and locations 
are determined via a T-expansion of the energy jfjj. We note that T-convergence of (1.1) 



to the functional ( |1.4[ ) with no screening and for fixed volume fractions was established by 
Ren and Wei in |43|, who also analyzed local minimizers of the sharp interface energy in 
the strong asymmetry regime in two space dimensions |42| . 

All these works are in the finite domain setting, while we are generally interested in 
the large volume (macroscopic) limit, i.e., the regime when the number of droplets tends 
to infinity. A rather comprehensive study of the behavior of the minimizers for the Ohta- 
Kawasaki energy in macroscopically large domains was recently performed in [36], still in 
the regime u close to —1. There the two-dimensional Ohta-Kawasaki energy was considered 
in the case when is a unit square with periodic boundary conditions. The interesting 
regime corresponds to the parameters e C 1 and 1 + u = 0(e 2//3 | In el 1 / 3 ) <C 1. It is shown 
in [36] that under these assumptions on the parameters and some technical assumptions 



on W, (1.4) gives the correct asymptotic limit of the minimal energy in (1.1). Moreover 



it is shown that when 5 := e 2 / 3 | lne| 1 / 3 (1 + u) becomes greater than a certain critical 



constant 5 C , the minimizers of E in (1.4) consist of 0(|lne|) simply connected, nearly 
round droplets of radius ~ S 1 / 3 ^ 1 / 3 1 lne| -1 / 3 , and uniformly distributed throughout the 
domain |36| . Thus, the following hierarchy of length scales is established in the considered 
regime: 

e < e 1/3 | lnel" 1 / 3 < | lne^ 1 / 2 < 1, (1.6) 

where the scales above correspond to the width of the interface, the radius of the droplets, 
the average distance between the droplets, and the screening length, respectively. The 



multiscale nature of the energy minimizing pattern is readily apparent from (1.6). 
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The analysis of 36 makes heavy use of the minimality condition for ( 1.4 ) and, in partic- 
ular, the Euler-Lagrange equation associated with the energy. One is thus naturally led to 
asking whether the qualitative properties of the minimizers established in [36] (roundness 
of the droplets, identical radii, uniform distribution) carry over to, e.g., almost minimizers 
of E, for which no Euler-Lagrange equation is available. More broadly, it is natural to ask 
how robust the properties of the energy minimizing patterns are with respect to various 
perturbations of the energy, for example, how the picture presented above is affected when 
the charge density u is spatially modulated. A natural way to approach these questions is 
via T-convergence. However, for a multiscale problem such as the one we are considering 
the proper setting for studying T-limits of the functionals in (1.1) or (1.4) is presently 
lacking. The purpose of this paper is to formulate such a setting and extract the leading 
order term in the T-expansion of the energy in ( 1 . 1 ) . In our forthcoming paper 24 , we 



obtain the next order term in the T-expansion, using the method of "lower bounds for 
2-scale energies" via T-convergence introduced in 
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For simplicity, as in |36| we consider the energy defined on a torus (a flat square with 
periodic boundary conditions). The main question for setting up the T-limit in the present 
context is to choose a suitable metric for T-convergence. In this paper we show that the right 
metric to consider is the weak convergence of measures for a suitably rescaled sequence of 
characteristic functions associated with droplets (see the next section for precise definitions 
and statements of theorems). Then, up to a rescaling, both the energy £ from (1.1) and 
E from (1.4) T-converge to a quadratic functional in terms of the limit measure, with the 
quadratic term generated by the screened Coulomb kernel from (1.5) and the linear term 
depending explicitly on 5 and k. 

We note that the obtained limit variational problem is strictly convex and its unique 
minimizer is a measure with constant density across the domain Q. In particular, this 
implies equidistribution of mass and energy for the minimizers of the diffuse interface 
energy 8 in (1.1) in the considered regime, which is a new result. In our companion paper 
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we further address the mutual arrangement of the droplets in the energy minimizing 
patterns, using the formalism developed recently for Ginzburg-Landau vortices [46] . In 
fact, the problem under consideration, and its mathematical treatment (here as well as 
in |24|), share several important features with the two-dimensional magnetic Ginzburg- 
Landau model from the theory of superconductivity 45 . There the role of droplets is 



played by the Ginzburg-Landau vortices, which in the appropriate limits also become 
uniformly distributed throughout the domain [441. We note, however, that the approach 



developed in [44 45 cannot be carried over directly to the problem under consideration, 
since the vortices are more rigid than their droplet counterparts: the topological degrees of 
the vortices are quantized and can only take integer values, while their counterparts here, 
which are the droplet volumes, are not. Thus we also have to consider the possibility of 
many very small droplets. Developing a control on the droplet volumes from above and 
below is one of the key ingredient of the proofs presented below, and relies on the control of 
their perimeter via the energy. Let us also mention yet another closely related system from 
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the study of ferromagnetism, where the role of vortices is played by the slender needle-like 
domains of opposite magnetization at the onset magnetization reversal |28| . 

Our paper is organized as follows. In Sec. [2j we introduce the considered scaling 
regime and state our main results; in Sec. [3] we prove the T-convergence result in the 
sharp interface setting; in Sec. [4] we prove the results on the characterization of almost 
minimizers of sharp interface energy; and in Sec. [5] we treat the T-limit for the case of the 
diffuse interface energy. 



Some notations. We use the notation (u £ ) G A to denote sequences of functions u £ G A 
as s = Eyi —7- 0, where A is an admissible class. We also use the notation fi G A4(£l) 
to denote a positive finite Radon measure dji on the domain tt. With a slight abuse of 
notation, we will often speak of [i as the "density" on Q. The symbols H 1 ^), BV{Q), 
C(Q) and i? _1 (fi) denote the usual Sobolev space, space of functions of bounded variation, 
space of continuous functions, and the dual of H 1 ^), respectively. 



2 Statement of results 

Throughout the rest of the paper the domain f2 is assumed to be a flat two-dimensional 
torus of side length £, i.e., when O = T? = [0,^) 2 , with periodic boundary conditions. 
Given the parameters k > 0, 5 > and £ > 0, for every e > we define 

u £ := -l + e^llnel 1 / 3 *. (2.1) 

Then the sharp interface version of the Ohta-Kawasaki energy (cf. ( |1.4| )) can be written as 

E £ [u} = - [ \Vu\dx+ 1 [ (u — n e )(— A + k 2 ) -1 (u — u £ ) dx, (2.2) 
2 JTj 2 h\ 

for all u G A, where 

A:=BV(Tj;{-l,l}). (2.3) 



We wish to understand the asymptotic properties of the energy E £ in (2.2) as e — > when 
all other parameters are fixed. We then relate our conclusions based on the study of this 
energy to its diffuse interface version, which under the same scaling assumptions takes the 
form 

S e [u] = j 2 \Vu\ 2 + W{u) + *(« - u £ ){-k)-\u - u £ )^j dx, (2.4) 



for all u G A £ , where 



A £ := ju G H l {Y}) : i J ^udx = if j . (2.5) 
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We note that at the level of the energy minimizers the relation between the two functionals 



was established in 36 



The energy E £ may alternatively be written in terms of the level sets of u. Indeed, 
the set VL + := {u = +1} is a set of finite perimeter (for precise definitions, see e.g. [3]), 
and when |<9f2 + | is sufficiently small, the set S7 + may be uniquely decomposed into an at 
most countable union of connected components flf, where the boundaries dQ.f of each 
connected component are Jordan curves which are essentially disjoint (after an extension 
to the whole of M? and up to negligible sets, see e.g. J2J Corollary 1 and Theorem 8]). In 
the context of T-convergence the sets Qf may be viewed as a suitable generalization of the 
droplets introduced earlier in the studies of energy minimizing patterns [36]. Note, however, 
that the sets 



36 



and may 



lack the regularity properties of the energy minimizers in 
in general be fairly ill-behaved (in particular, they do not have to be simply connected). 
Nevertheless, they are fundamental for the description of the low energy states associated 
with E £ and, in particular, will be shown to be close, in some average sense, to disks of 
prescribed radii for almost minimizers of energy. 

In terms of the droplets, from the above discussion we have 



u 



■1 



2 E^ 



(2.6) 



where Xo+ are the characteristic functions of 0^~. Inserting this into (2.2), expressing the 
result via G that solves 

-AG(x) + k 2 G(x) = S(x) in Tj, (2.7) 

expanding all the terms and using the fact that J" t2 G(x)dx = k~ 2 , we arrive at (see 
also [36] ) 

l 2 {l + u e ) 2 



E £ [u] 



2k 2 



+ J2{ £ \ dn t\ -2k~ 2 (1 + u £ )|S!+|} +2^ / f + G{x-y)dxdy, (2.8) 

where we took into account the translational symmetry of the problem in T». Moreover, 
since the optimal configurations for are expected to consist of droplets of size of order 
e 1//3 | lne|~ 1//3 (see (1.6) and the discussion around), it is convenient to introduce the rescaled 



area and perimeter of each droplet: 

Ai := e~ 2/3 \ lne\ 2/3 \nf\, Pi := e~ 1/3 \ lne\ 1/3 \dnf\. (2.9) 
Similarly, let us introduce the suitably rescaled measure d\x associated with the droplets: 

dfi(x) :=e- 2/3 \lne\- 1/3 J2x n +(x)dx= 1 £" 2/3 | lne\- 1/3 (l + u) dx. (2.10) 

i 2 
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Then the energy E £ [u] may be rewritten as 



E £ [ U ] = e^\lne\ 2 / 3 (^ + E £ [ U ] ] j, 



(2.11) 



where 



E £ \u] 



2 ^A! ) + 2 



G{x-y)dn{x)dn{y). (2.12) 



Observe that for the nontrivial minimizers we know from 36 that E e = 0(1), A4 = O(l) 
and Pi = O(l) (and even more precisely Ai ~ 3 2//3 7r and Pi ~ 2 • 3 1//3 7r), the number of 
droplets is N = 0(|lne|), and \i closely approximates the sum of Dirac masses at the 
droplet centers with weights of order | lne| _1 . Observe also that if (u £ ) S A is such that 
lim sup e ^ E £ [u £ ] < +00, we have 



lim sup — r Pf < +00, lim sup — r A £ < +00, 

E ^n \ms\ *— J r^n In e ^ 



(2.13) 



where {^4f}, {Pf} and fi £ are associated with u £ . Indeed, by (2.9) and (2.10) we have 



(2.14) 



and thus, from (2.12) we obtain 
E £ [u] > -% 



d[i £ + 2 



> 



26 

Jt\ 



G(x -y)dfi £ (x)dfi £ {y) 
/ \ 2 



d/j £ + 



K 



2 £ 2 



d[i £ 



(2.15) 



where we used the well-known fact that for any u £ L (Tj) we have 



G(x — y)u(x)u(y)dx dy > 



i 2 



G(x) dx I I / u(x) dx 



which is readily obtained by Fourier transform. Then it follows that 

1/2 



T 2 



dfi £ < 5t + k£ 



max(0,E £ [u £ }) 



(2.16) 



(2.17) 



implying both statements in (2.13). 



S 



We now state our T-convergence result, which is obtained for configurations (u £ ) that 
obey the optimal energy scaling, i.e. when E"^?/] remains bounded as e — > 0. The result 
is obtained with the help of the framework established in 44 , where an analogous result 
for the Ginzburg-Landau functional of superconductivity was obtained. What we show 
is that the limit functional E° depends only on the limit density \x of the droplets (more 
precisely, on a general positive measure d\i). In passing to the limit the second term in 
(2.12) remains unchanged, while the first term is converted into a term proportional to 
the integral of the measure. The proportionality constant is non-trivially determined by 
the optimal droplet profile that will be discussed later on. We give the statement of the 
result in terms of the original screened sharp interface energy E £ , which is defined in terms 
of u G A. In the proof, we work instead with the equivalent energy E £ , which is defined 
through {Af}, {if} and fi £ (cf. pll] ) and fl2l5| |). 

Theorem 1. (r-convergence of E £ ) Fix 6 > 0, k > and £ > 0, and let E e be defined 
by (2.2) with u £ given by (2.1). Then, as e — > we have that 

e-^\ lne\-^Et 4 E° \p] := ^ + ( 3 2 / 3 - %) [ d„ + 2 / / G(x - y)d^x)d„(y), 

where fi G Ai(Tj) n H^ 1 {Tf). More precisely, we have 
i) (Lower Bound) Let (u £ ) £ A be such that 

limsupe- 4/3 | lne\- 2/3 E £ [u £ ] < +oo, (2.18) 

e^O 

let 

dfi £ (x) := ^e~ 2/3 | lne|" 1/3 (l + u £ (x))dx, (2.19) 

and let v £ satisfy 

- Av £ + k 2 v £ = // in Tj. (2.20) 
Then, up to extraction of a subsequence, we have 

ff -± ii in {C(Tj)Y, v £ -^v in H l (Tf), 
ase^O, where fj, G M{Tf) n H- l (Tf) and v G fl^Tf) satisfy 

-Av + k 2 v = h in T 2 e . (2.21) 

Moreover, we have 

liminfe- 4 / 3 |lner 2 / 3 E £ [n £ ] > E°\u]. 
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ii) (Upper Bound) Conversely, given fj, G ,M(T 2 ) n H 1 (T|) and v G H (Tj) solving 



2.21), there exist (u £ ) G *4 suc/t that for the corresponding [i £ , v £ as in (2.19) and 



(2.20) we have 



as e — > 0, and 



if-LfJkin (C(Tf))*, /-i;in if^Tf), 

£^0 



We note that the limit energy E° obtained in Theorem [T] may be viewed as the ho- 
mogenized (or mean field) version of the non-local part of the energy in the definition of 
E £ associated with the limit charge density [i of the droplets, plus a term associated with 
the self-energy of the droplets. The functional E° is strictly convex, so there exists a 
unique minimizer p, G .A/f (T|) n H~ (T|) of E°, which is easily seen to be either fx = 

\(8- 53 2/3 k 2 ) otherwise. The latter can also be seen immediately 
below, which gives a local characterization of the limit energy E°. 



for 5 < |3 2 / 3 k 2 or fl 
from Remark 



2.1 



Remark 2.1. The potential v from ( 2.21[ ) can also be expressed as v(x) = J" t2 G(c 



y)d[i(y). The limit energy E° in Theorem^ becomes local when written in terms of that 



limit potential v: 



2e2 



5H 
2k/ 



+ 



3 2 /V 



28 



v dx + 2 



Vv\ 



K, 2 v 2 ^jdx. 



(2.22) 



Also, by the usual properties of T-convergence IB], the optimal density \i above is 
exhibited by the minimizers of E £ in the limit e — > 0, in agreement with 36, Theorem 2.2]: 



Corollary 2.2. For given 5 > 0, K > and £ > 0, let u £ be given by (2.1 ) and let (u £ ) G A 



be minimizers of E £ . Then, if /J, 6 is given by (2.19), as e — > we have 



e- 4/3 |lner 2/3 r 2 min£; £ 




7?(2<$-<y, 



(2.23) 



w/ien 5 < <5 C or 5 > <5 C , respectively, with 5 C 



In particular, since the minimal energy scales with the area of T», it is an extensive quantity. 

We next give the definition of almost minimizers with prescribed limit density, for 
which a number of further results may be obtained. These can be viewed, e.g., as almost 
minimizers of E £ in the presence of an external potential. We note that in view of strict 
convexity of E°, minimizing E°[fi] — J t2 ip(x)d[i(x) for a given ip G ^(ir 2 ) one obtains 
a one-to-one correspondence between the minimizing density li and the potential (p. It 
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then makes sense to talk about almost minimizers of the energy E £ with prescribed limit 
density [i by viewing them as almost minimizers of E £ — Jj 2 (pd/j, £ . Also, observe that 



almost minimizers with the particular prescribed density \i from Corollary 2.2 are simply 
almost minimizers of E £ . Below we give a precise definition. 

Definition 2.3. For given 5 > 0, k > 0, £ > and a given ^ G M(Tj) n H^ijf), we 
will call every recovery sequence (u £ ) G A in Theorem^ii) almost minimizers of E e with 
prescribed limit density [i. 

For almost minimizers with prescribed limit density, we show that in the limit e — > 
most of the droplets, with the exception of possibly many tiny droplets comprising a van- 
ishing fraction of the total droplet area, converge to disks of radius r = 3 1//3 e 1//3 | lne| -1 / 3 . 
More precisely, we have the following result. 

Theorem 2. Let (u £ ) G A be a sequence of almost minimizers of E £ with prescribed limit 
density \i. For every 7 G (0, 1) define the set I £ := {i G N : 3 2 / 3 ttj < A\ < 3 2 / 3 7T7 -1 }. 
Then 

lS5|hb|E(^-VS3f)=0, (2.24) 

i 

>i?„iib|E(^-^) 2 =°' <" 5 > 



lim M = 0, 2.26) 

e^o lne 

1 1 im 



where {A £ } and {P £ } are given by (2.9) with u 



Note that we may use the isoperimetric deficit terms present in (2.24) to control the 
Fraenkel asymmetry of the droplets. The Fraenkel asymmetry measures the deviation of 
the set E from the ball of the same area that best approximates E and is defined for any 
Borel set E C M 2 by 

/ s \EAB\ 
a(E) = min 1 1 , (2.27) 
\E\ 

where the minimum is taken over all balls B C M 2 with \B\ = \E\, and A denotes the 
symmetric difference between sets. Note that the following sharp quantitative isoperimetric 
inequality holds for a(E) [2l] : 

per(E) - s/krr\E\ > Ca(E)^/\E~\, (2.28) 



with some universal constant C > 0. As a direct consequence of Theorem [2j and ( |2.14| ) , we 
then have the following result. 
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Corollary 2.4. Under the assumptions of Theorem^ when Jj 2 dfj, > we have 



3 2 /3 7r |je| r -, 

lim - , 7 = / du, lim — > a(H+) = 0, (2.29) 

II •> il £ 171 jg/s 



where \I £ \ denotes the cardinality of I £ . 

This result generalizes the one in [36] , where it was found that in the case of the minimizers 
a// the droplets are uniformly close to disks of the optimal radius r = 3 1//3 e 1//3 | lne| -1 / 3 . 
What we showed here is that this result holds for almost all droplets in the case of almost 
minimizers, in the sense that in the limit almost all the mass concentrates in the droplets 
of optimal area and vanishing isoperimetric deficit. We note that the density \i is also 
the limit of the number density of the droplets, up to a normalization constant, once the 
droplets of vanishing area have been discarded. 

The result that almost all droplets in almost minimizers with prescribed limit density 
have asymptotically the same size, even if the limit density is not constant in T| appears 
to be quite surprising. In fact, this observation is consistent with the expectation that 
quantum mechanical charged particle systems form Wigner crystals at low particle den- 
sities 25 32 54 . Let us point out that the Ohta-Kawasaki energy £ £ bares a number of 



similarities with the classical Thomas-Fermi-Dirac-Von Weizsacker model arising in the 



context of density functional theory of quantum systems (see e.g. 30-32]). 

We now turn to relating the results obtained so far for the screened sharp interface 
energy E £ to the original diffuse interface energy £ £ . On the level of the minimal energy, 
the asymptotic equivalence of the energies in the considered regime, namely, that for every 
5 > 

(1 - 5) m\nE £ < mm£ £ <(!-<$) mmE e (2.30) 



for e -C 1 was established in |36[ Theorem 2.3]. The main idea of the proof in [36] is 
for a given function u e S A £ to establish an approximate lower bound for £ £ [ii £ ] in terms 
of (1 — 5)E £ [u £ ] for some u e £ A, with 5 > which can be chosen arbitrarily small for 
£< 1. The matching approximate upper bound is then obtained by a suitable lifting of 
the minimizer u £ G A of E £ into A £ . 

Here we show that the procedure outlined above may also be applied to almost mini- 



mizers of E in a suitably modified version of Definition 2.3 involving £ £ . We note right 
away, however, that it is not possible to simply replace E £ with £ £ in Definition 2.3 The 
reason for this is the presence of the mass constraint in the definition of the admissible class 
A 6 for £ £ . This implies, for example, that any sequence of almost minimizers (u £ ) £ A £ 
of £ £ must satisfy l~ 2 J* T 2 dfi £ = ^5, while, according to Corollary 2.2, for sequences of 

almost minimizers (u £ ) G A of E £ we have £~ 2 J t2 d[i £ — >• p ^ ^5. This phenomenon is 
intimately related to the effect of screening of the Coulombic field from the droplets by the 
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u 

1.0 : 
0.5 - 



-0.5 - 



Figure 2: A qualitative form of the u-profile for a single droplet from the Euler-Lagrange 
equation associated with E. The horizontal line shows the level corresponding to u. Charge 
is transferred from the region where u < u (depletion shown in green) to the region where 
u > u (excess shown in orange). 



compensating charges that move into their vicinity [35] . For a single radially symmetric 
droplet the solution of the Euler-Lagrange equation associated with E £ will have the form 
shown in Fig. [2j 

In order to be able to extract the limit behavior of the energy, we need to take into 
consideration the redistribution of charge discussed above and define almost minimizers 
with prescribed limit density that belong to A £ and for which the screening charges are re- 
moved from the consideration of convergence to the limit density. Hence, given a candidate 
function u £ 6 A 6 , we define a new function 

e ( \ /+!> u £ (x)>0, 
Un(x) := < (2.31) 
\-l, u £ (x) <0, V ; 

whose jump set coincides with the zero level set of u £ . This introduces a nonlinear filter- 
ing operation that eliminates the effect of the small deviations of u £ from ±1 in almost 
minimizers on the limit density (compare also with [28]). The density is now defined as 

dji% := \e- 2/z \ lne\- 1/3 {l + u £ (x))dx. (2.32) 

We can follow the ideas of 36 to establish an analog of Theorem [T] for the diffuse 
interface energy. To avoid many technical assumptions, we formulate the result for a 
specific choice of W. A general result may easily be reconstructed. Also, we make a 
technical assumption to avoid dealing with the case limsup e _>. H^H^occ^) > 1, when spiky 
configurations in which \u £ \ significantly exceeds 1 in regions of vanishing size may appear. 



We note that this condition is satisfied by the minimizers of £ £ 36, Proposition 4.1] 
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Theorem 3. (r-convergence of £ £ ) Fix 5 > and £ > 0, and let £ £ be defined by (2.4) 
with W(u) = ^(l — u 2 ) 2 and u £ given by (2.1). Then, as e —¥ we have that 



2 1)2 



2k 2 



3 2/3 



2(5 



dfi + 2 



G{x - y)dfj,(x)dfj,(y), 



where fi G .A/f (T 2 ) n -ff 1 (T^) and k = |. More precisely, we have 



i) (Lower Bound) Let (u £ ) G >4 e 6e suc/i i/icrf limsup e _ s>0 ll ue |lz,°°(Tf) — 1 an d 

limsup£~ 4/3 |ln{r|~ 2/3 £ £ [?/] < +oo, 

£^0 



(2.33) 



and let Hq(x) be defined by (2.31) and (2.32). 
Then, up to extraction of subsequences, we have 

as e — > 0, where fi G M.(Y 2 )C\H~ l (Y 2 ) . Moreover, we have limsup^Q II u£ |Il°°(t|) = 1 
and 

liminfe- 4 / 3 |lner 2 / 3 £: £ [u £ ] > E°U. 

£->0 

ii) (Upper Bound) Conversely, given /i G .A/f (T 2 ) n H~ l (T 2 ), there exist (u £ ) G A £ such 
that limsup^^Q H'u 5 ||z,°°(T|) = 1 and for fi^ defined by ( 2.31| ) and (2.32) we have 



Hl^vin (C(T 2 ))*, 



as e — > 0, and 



limsup^ 4/3 |lne|- 2 / 3 ^[u e ] < E° [fj] . 



Based on the result of Theorem [3| we have the following analog of Corollary 2.2 for the 
diffuse interface energy £ £ . 



Corollary 2.5. For given 5 > and t > 0, let u £ be given by (2.1) and let (u £ ) G A £ be 
minimizers of £ £ defined in (2.4) with W(u) = J>(1 — u 2 ) 2 . Then, if Uq and are defined 



via (2.31) and (2.32), respectively, as e — )• we have 



Mo 







in 



when 5 < 5 C or 5 > 5 C , respectively, with 5 C : 



(C(T 2 ))*, e - 4 / 3 |lne|- 2 / 3 r 2 min^ -+ { 2 f 



%{25 -S c 



(2.34) 



§3 2 / 3 k 2 and k = §. 



Remark 2.6. T/ie choice of the zero level set ofu £ in the definition of the truncated version 
Uq of u £ in ( |2.31 ) was arbitrary. We could equivalently use the level set {u £ = c} for any 
c G (—1, 1) fixed. 
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3 Proof of Theorem [T] 



Throughout all the proofs below, the values of A\ and Pf are always the rescaled areas 
and perimeters, defined in ( |2.9[ ), of the connected components flf of il + = {u = +1} for a 
given u = u £ , as in Sec. [2| The presentation is clarified by working with the rescaled energy 
E £ defined by (2.12) rather than E 6 directly. We begin by proving Part i) of Theorem [TJ 



the lower bound. 

3.1 Proof of lower bound, Theorem [I] i) 

Step 1 : Estimate of E £ in terms of AS and Pf . 

First, for a fixed 7 £ (0,1) we define a truncated rescaled droplet area: 
; £ . = Uf, if At < S 2 / 3 ^- 1 

*' \(3 2 /3 W7 -l)l/2|^||l/2 if A e > 32/3^-^ I- ■ J 

and the isoperimetric deficit 

^ :=n^E - v^t ) > 0, (3.2) 

which will be used throughout the proof. The purpose of defining the truncated droplet 



area in (3.1) will become clear later. We will make repeated use of the basic fact that the 
diameter of an indecomposable set of finite perimeter in the plane is essentially controlled 
by its perimeter (i.e., modulo a set of measure zero): 

essdiam u> < *|<9w| Vu Cl 2 , < |w| < 00. (3.3) 

The proof follows, e.g., from |2j Theorem 7 and Lemma 4] by noting that in view of [2j 
Proposition 6(ii)] it is sufficient to consider only simple sets [2j Definition 3]. 
We start by writing // = £V /i? , with 

dnl{x) := e" 2/3 | lne|- 1 / 3 Xo+(2;)^, (3-4) 

i 

where are the connected components of S7 + = {u e = +1}, and the index e was omit- 
ted from to avoid cumbersome notation. For small enough e this is justified by the 
discussion of Sec. [2j in view of the fact that for some C > we have 

\dn + \ < e^E^u 6 ] < Ce 1 / 3 ! lne| 2 / 3 , (3.5) 



15 



so \dQf\ <C £ whenever e <C 1. In particular, we can apply (3.3) to W~ on Tj to obtain 



essdiam flf < -\dO,f\, 



(3.6) 



when e is sufficiently small. 

For a fixed p > we introduce the "far field truncation" G p G C°°(T^) of the Green's 
function G: 



G p (x-y) :=G{x-y)4> p {\x-y\) V(x, y) G Tj x Tf, 



(3.7) 



where t/>p G C°°(R) is a monotonically increasing cutoff function such that 4> p (t) = for 
all £ < 7^p and ct) p {t) = 1 for all t > p. Then, for sufficiently small e we have |<9f2^| < p in 



view of (3.5), and from (2.12) we obtain 



+ 2 E // G{x-y)d l x\{x)d l 4{y) 
+ 2 [[ G p {x-y)dp e {x)dpi £ {y), 



(3.8) 



where we used (3.6) and the positivity of G (cf. e.g. |36|), and here and everywhere below 
we omit T 2 , x T|. as the domain of integration for double integrals to simplify the notation. 
We recall that the Green's function for —A + k? on T| can be written as G(x — y) = 



— ^ In \x — y\ + 0(\x — y\) [36| . With the help of this fact, together with ( |3.5| ) and (3.6) 
for e sufficiently small we have the following estimate for the self-interaction energy: 



^sclf 



i 

= -^Lf2 E f u+ 1 + ( ln (k 1/3 i ln£ i 2/3 i- -y\) + c ) 



for some C > independent of e, where in equation (3.9) we have rescaled coordinates 
x = e _1//3 | lne| 1//3 x, y = e _1 / 3 | lnel 1 / 3 and introduced the rescaled versions fij of 
Expanding the logarithm in (3.9) and using (2.13) and (3.6), we obtain that E^ ei{ can be 
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bounded from below as follows: 



-^self > 



> 



> 



i 

1 « 



3tt 



c 



In | lne| 
I me l 

In | lne| 
I me l 

In | lne| 
I In el 



1 



In \x — y\dx( 



1 



7r| lnel 



InH 



(3.10) 



for some C > independent of e (which changes from line to line) . 



Now observe that the term in parentheses appearing in the right-hand side of (3.10) 



is positive for e sufficiently small. Using this and the fact that A\ > Af, from (3.10) we 
obtain 



— Vli £ l 2 f- 



C 



In | lne| 
I lnel 



(3.11) 



where C > is a constant independent of e. It is also clear from the definition of A\ that 
there exists a constant c > such that 



if I 2 < cAS. 



(3.12) 



Combining this inequality with (3.11) and choosing any 5 > 0, for e small enough we have 



5 > Cc ^ngp and, therefore, from (3.8) we obtain 



+ 2 if G p (x-y)dn e (xW £ (y). 



(3.13) 



Step 2: Optimization over A 6 , . 



Focusing on the second term in the right-hand side of (3.13), we define 

f(x) := h/* + Jl Xj (3.14) 

and observe that / is strictly convex and attains its minimum of 3 2 / 3 at x = 3 2 / 3 tt, with 

3a/tt 



/"(*) 



2x 5/2 ' 



(3.15) 
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We claim that we can bound the second term in the right-hand side of (3.13) from below 
by the sum I + II + III of the following three terms: 



1 



II 



III 



1 



3 2/3_ 
7 5/2 



7 i ' Af>3 2 / 3 7r7-! 



1)A: 



lne|4vr 2 -3 2 / 3 



£ Af(Af-3 2 / 3 vr) 5 



Af<3 2 / 3 7T7 



1 7 7 / 2 

In el 47r 



E 



(Af-3 2 / 3 vr) 2 - 



(3.16) 
(3.17) 
(3.18) 



3 2 / 3 7r 7 <ylf<3 2 / 3 7r7- 1 

Before proving this, observe that defining 



"> e i-n^( 32/3 -i?-«)s>*-'// 



GJx-y)dif{x)dif{y), (3.19) 



we have from ( 3.13[ ) and (3.16)-( |3.18[ ) that if 1^ is as in Theorem[2j then 



M£ ^ HnV E ^ + riti E(^f - 32/3 -) 2 + 3* ^ v ^ (o, } 



3^' 



(3.20) 



for some constants c\, C2 > depending only on 7. 

We now argue in favor of the lower bound based on (3.16)-(3.18). First observe that 
by dO) we have for all Af > tf^n-y' 1 : 



1 



25 



+ 5) > (V 3 - ^ - 6) A\ + 3 2 / 3 (3- 1 7 _1 - l)Af . (3.21) 



When Af < 3 2//3 7T7 1 , which corresponds to both (3.17) and (3.18), we utilize convexity of 
/ and pOo): 



25 



+ 8\JE = M 



2^fn 1 25 / 

— 7= H n - 5 

37T fc 



Af /(Af 



25 



(3.22) 



> (3 2/3 - ^ - Af + \A\f (S 2 / 3 ^- 1 ) (Af - 3 2 / 3 vr) 2 , 



where the last line follows from the second order Taylor formula for f(x) about x = 3 2 / 3 7r 
and the fact that f"(x) is decreasing. Combining (3.19), (3.21) and (3.22) yields M £ > 
I + 11 + III. 
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Now using (3.19) and (3.20) with 7 sufficiently small, we deduce that 

, p (x - y)dfj, £ (x)d/j, £ (y). 



^«*]^( 32/3 -^H*>' +2 // G ' 



(3.23) 



Step 3: Passage to the limit. 



We may in fact conclude from (2.12), (2.13), (2.18) and (2.20), that 



limsup / \Vv £ \ 2 + k 2 \v £ \ 2 < +00, 



while ([i £ ) is bounded in the sense of measures from (2.13) and (2.14). 
Consequently, up to a subsequence 

v £ v in H 1 ^), 
// ^ in C(Tf), 

where 

— Av + k 2 v = fl 



(3.24) 



(3.25) 
(3.26) 

(3.27) 



holds in the distributional sense. Now passing to the limit in (3.23) and recalling (2.14) 
we obtain 



liminf E £ [u £ ] > ( 3 2/3 



25 



dfi + 2 // G p (x-y)dfi(x)dfi(y), 



(3.28) 



using continuity of G p . An application of the monotone convergence theorem then yields 



liminf E £ [u £ ] > U 2 > z - ^ J / <ln + 2 / / G(:r - //)r///(.r )<///(//). 



upon sending p — >• and then 5—^0. 



(3.29) 
□ 



We now argue in favor of the corresponding upper bound in Theorem [T] The construc- 
tion resembles quite closely that of the vortex construction in 44 for the two dimensional 



Ginzburg-Landau functional and indeed we borrow several ideas from that proof and oc- 
casionally refer the reader to that paper for details. 
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3.2 Proof of the Upper Bound, Theorem [T] ii) 



As in the proof of the lower bound, we set dfj,f(x) as in ( |3.4[ ), so that \x £ = ^ \i\ . Using the 
approximation argument of Proposition II. 2 in [44] we may assume that the limit measure 
fi satisfies 



dfj,(x) = g(x)dx, c < g < C, 



(3.30) 



for some C > c > 0. 



Step 1: Construction of the configuration. 

We claim that for e sufficiently small it is possible to place a total of N(e) disjoint spherical 
droplets, where 



N(e) 



1 |lne| 

32/3^" 



MT, 2 ) + o(|ln £ |), 



with centers {aj} in T| and radius 



and satisfying for all i 7^ j 



d(e) := min |oj — dj\ > 



(3.31) 



(3.32) 



(3.33) 



for some constant C > depending only on /j,. Indeed, given \x satisfying (3.30), for e 
sufficiently small we can partition T» into disjoint squares {K{\ of side length 5 e > 
(hereafter simply denoted 5) satisfying 



lner 1 / 2 « < y< 1. 



In each K{ we place 



1 |lne| 
32/3 vr 



(3.34) 



(3.35) 



points dj (here m = \x\ denotes the smallest integer m < x) satisfying (3.33) and in 
addition 

dist fadKi) > - 7 |=, N(e) :=^N Ki . (3.36) 

As argued in 44 j, our ability to do this follows from the estimate: 

c5 2 < n(Ki) < C5 2 , 



(3.37) 
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which follows from (3.30) together with (3.34). We finally define our configuration u £ by 



setting the connected components £lf of f2 + = {u 6 = +1} to be balls of radius r from 
(3.32) centered at aj, i.e. := B r (a{). We set u £ = — 1 in the complement of these balls. 
With these choices we have 



- E[u£] = 2vr • S^Nje) _ 2^.^N{e)5 + , y _ 



In el 



lnel/t 2 



G(x - y)dp e {x)dp e {y) 



^p(T 2 e ) - ^(T 2 ) + 2 J! G{x - y)d l f{x)d l f{y) + o(l). (3.38) 



The main point of the rest of the proof is to show that the integral term in (3.38 ) converges 



to J J G{x — y)dp{x)dp(y) + 3 -1 / 3 J dp, with the non-trivial last term coming from the self- 
interaction of the droplets. To prove that these are the only contributions to the limit 
energy, we need to use the fact that the droplets cannot concentrate too much as e — > 0. 

Step 2: Convergence of the configurations. 

Defining p e as before, it is clear from the construction that 

// - // in (C(T 2 ))*. (3.39) 
Fix p > sufficiently small (depending only on n and £) and consider G p {x — y) defined as 



in (3.7). By continuity of G p in Tj we have 



lim 

e-»o J J 



/ / G p( x ~ y)dp £ (x)dp e (y) = G p (x - y)dp(x)dp(y). 



(3.40) 



Now, let I p be the collection of indices such that < |o, — Oj| < p. Then for small 
enough e we can write 



(G(x -y)- G p {x - y))dp £ {x)dp £ (y) 

m ,, ,, 
<J2 G(x- y)dpKx)dp!(y) + £ G(x - y)dp!(x)dp £ (y) (3.41) 
i=i 7 7 (i,j)ei P 

N(e) 

I 

< 



67r| In el ^— ' 1 11 ' I In el ' |lne| 2 



for some C,C' > independent of e or p, where A\ = 3 2 / 3 7r and we expanded the Green's 
function as in (|3.9[) in the proof of the lower bound. Now, for k = 1,2,..., K p (e), with 



K p (e) := [p/d(e)\, let 7™ C I p be disjoint sets consisting of all indices such that 

kd(e) < \ai — a,j\ < (k + l)cZ(e). Since by the result on optimal packing density of disks in 
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the plane |19| we have |J^| < ckN(e) for some universal c > (here again denotes the 



cardinality of 1%), in view of (3.31) it holds that 



lne 



K P {e) 



^ £ a\a £ indist (n+,n+) < £ 



fc=l 



< 



2CiV(e) 
lne| 2 (i 2 (e) 



" t|lnt|dt<C , r |h ! d(£)l +P 2 \lnp\ ) <2CV|lnp|, 



lne 



(3.42) 



for some C, C > independent of e or p, when e and p are sufficiently small. Therefore, 



from (3.32) and (3.41) we obtain 



limsup / / (G(x -y)- G p (x - y))dp £ (x)dp £ (y) < 2" 1 • 3" 1/3 + o(p). (3.43) 



Finally combining (3.43) with (3.38) and (3.40), upon sending e — > 0, then p — > and 



applying the monotone convergence theorem we have 

lim^ £ [n £ ] < (V /3 -^J J dp + 2 JJ G(x- y)dp(x)dp(y), 



(3.44) 



as required. The fact that v e v follows from (3.39) and the uniform bounds just demon- 



strated on the terms involving the Green's function in (3.38), from which it follows that 
(2.21) is satisfied distributionally. □ 



4 Proof of Theorem 



In the proof of Sec. [5J we have in fact established Theorem [2j which is clear by (3.20). 
Indeed, we have for a sequence of almost minimizers (u £ ): 



]im E £ [u £ ] -E \u] = 0. 

£->0 



(4.1) 



Observing that M £ defined in (3.19) does not contribute to Eo[p], we have e stabli shed that 
M £ — > as e — > for any 7 < 3 and, as a consequence, we obtain ( 2.24 )-( 2.26) for, say, 
7 = g. Then it is easy to see from the definition of I £ that the statement of the Theorem, 
in fact, holds for any 7 £ (0,1). □ 



5 Proof of Theorem [3] 

We now turn to the proof of Theorem [3] extending the result of Theorem [T] for the sharp 
interface energy E £ to the diffuse interface energy E £ . The proof proceeds by a refinement 
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of the ideas of 36 , Sec. 4] to establish matching upper and lower bounds for £ £ in terms 
of E £ for sequences with bounded energy. 

Step 1: Approximate lower bound. 



In the following, it is convenient to rewrite the energy (2.4) in an equivalen form 

.2 



£ £ [u e 



%\Vu £ \ 2 + W(u £ ) + l\Vv £ \ 2 
2 \ 2 2 

e 



dx, 



-Av £ 



u 



u £ , 



v £ dx = 0. 



(5.1) 



Fix any 5 £ (0, 1) and consider a sequence (u £ ) £ A £ such that limsup e _> ||^ £ ||z°°(t 2 ) — 1 
and £ £ [u £ ] < Ce^ 3 \ lne| 2 / 3 for some C > independent of e. Then we claim that 



limsup ||u e ||^ 



°°(Tf) 



lim ||t; £ II j oo/ T 2^ 



0. 



Indeed, for the first statement we have from the definition of £ £ in (2.4) that 
|^|<Ce 4 / 3 |lne| 2 / 3 <5- 2 , 



O 5 - 



{-l + <5 < u £ < 1-5}, 



(5.2) 



(5.3) 



for some C > independent of e. Hence, in particular, limsup^Q ||^ e ||ioo(x 2 ) — 1> proving 
the first statement of (5.2). To prove the second statement in (5.2), we note that by 
standard elliptic theory 22 we have ||w e ||vi/ 2 ^(T 2 ) — C f° r an y P > 2 and some C" > 
independent of e and, hence, by Sobolev imbedding ||V-u £ || L oo( T 2) < C" for some C" > 
independent of e as well. Therefore, applying Poincare's inequality, we obtain 



Ce 4/3 |lne| 2 / 3 > £ £ [u £ ) > C f 



v £ \ 2 dx > C"\\v 



"IL.e||4 



Il°°(t|)' 



for some C , C" > independent of e, yielding the claim. 



In view of (5.2), for small enough e we have ||it e ||nx>(x 2 ) — 1 + ^ 3 an d \\ vS 



lz°°(T 2 ) 



(5.4) 



<<5 3 , 



and by the assumption on energy we may further assume that £ £ [u £ ] < 5 12 . Therefore, 
by 36, Proposition 4.2] there exists a function Uq E A such that 

£ £ [u £ ]>(l-5 1/2 )E £ [u £ }. (5.5) 

In particular, (uq) satisfy the assumptions of Theorem [TJ and, therefore, upon extraction 
of subsequences we have jip/ie M(Tf) n H~ l (Y 2 g ) in (C(T 2 ))*, where 



d/xg(s) := ±e~ 2/3 |m£|-^ 3 (l + ^(x))cfa. 



(5.6) 



Furthermore, recalling that by construction the jump set of Uq coincides with the level set 
{u £ = c} for some c S (—1 + 6,1 — 5), see the proof of [36[ Lemma 4.1], from (5.3) we have 



l^-^0llL 1( T?)<^ 4/3 | lne | 2/3r2 ' ( 5 - 7 ) 
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where Uq is given by ( |2.31 ), for some C > independent of e (recall that the jump set of 
Uq is the zero level set {u £ = 0}). Comparing ( |5.7[ ) with (5.6), we then see that Hq — fj, in 
(C(T 2 ))* as well. The result of part (i) of Theorem [3] then follows by arbitrariness of 5 > 
via a diagonal process. □ 

Step 2: Approximate upper bound. 

We now construct the approximate upper bounds for a suitable lifting of the recovery 
sequences in the proof of Theorem [jjii) to A £ . Let (mq) E A be a recovery sequence 
constructed in Sec. 3.2 This sequence consists of circular droplets of the optimal radius 
r = S 1 / 3 ^ 1 / 3 ! lne| -1 / 3 > e 1 / 2 and mutual distance d > C\ lne| -1 / 2 > e 1 / 2 , for some C > 
independent of e. In addition, since 



£ e [it§] = -/ \Vu £ \dx + - / |V£T + k>T dx < C£ 4/3 |lne|^ 3 , (5.8) 
where £ e (x) = | l2 G(i — y){u £ Q{y) — u £ )dy, for some C > independent of e, by the 



argument of (5.4) one can see that lim e _>o ||^ £ ||.l°°(t 2 ) = 0- Therefore, for any 5 G (0,1) 
and e > sufficiently small we have ||£ £ ||l°o(t2) < <5 and .E e [uQ] — <5 5 ^ 2 - So we can 
apply [36j Proposition 4.3] to obtain a function u £ G .A e such that 

£ e [u £ }<(l + 5 l ' 2 )E e [ul]. (5.9) 

Furthermore, by the construction of u £ (see |36[ Eqs. (4.31)-(4.33)]) and arbitrariness of 
5 > 0, we also have limsup 



£->0 



lL°°(Tf) 



1, and 



< Ce 4/3 |lne| 



2/3 



(5.10) 



for some C > independent of e, where tig is given by (2.31), and we used (5.8). Hence 
^0 A* 

arbitrariness of 5 > via a diagonal process. 



lim e _>o m (C(^))*- The resu lt of part (ii) of Theorem [3] again follows by 

□ 

Remark 5.1. It is possible to chose 5 = e a for a > sufficiently small in the arguments 
of the proof of Theorem^ Therefore, given a sequence of almost minimizers (u £ ) £ A £ of 
£ £ and the corresponding sequence (tig) £ A of almost minimizers of E £ , one has 



£" 4/3 | ln£\-^ 3 £ £ [u £ ] = £~ 4/3 | lne\- 2/3 E £ [u £ ] + 0(e c 
for some a <C 1, as e — >• 0. 



(5.11) 
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